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A phenomenon analogous to the conical refraction widely known in the crystalooptics and crys- 
taloacoustics is discovered for the magnetohydrodynamical waves in the coUisionless plasma with 
anisotropic thermal pressure. Angle of the conical refraction is calculated for the medium under 
study which is predicted to be 18° 26'. Possible experimental corroborating of the discovered phe- 
nomenon is discussed. 
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> ' I. INTRODUCTION 

o 

, A phenomenon of the conical refraction — a pecuhar kind of hght refraction in biaxial crystals taking place when a 
00 ■ direction of a light beam coincides with some of the optical axes (binormals) of the crystal, was predicted theoretically 
' by W.R. Hamilton even in 1832. He discovered this phenomenon basing on the Huygens-Fresnel principle (formulated 
17 years earlier) applied for the biaxial crystals. In 1833 H. Lloyd proved experimentally existence of the phenomenon. 
" ' ' Observation of the inner conical refraction in biaxial crystal is possible if the parallel, narrow beam of natural (non- 
, polarized) light falls along the optical axis (binormal) on a plate cut out from the crystal perpendicularly to the 
binormal. Propagating through the crystal the beam will form an empty cone with permanently changing linear 
, polarization. Further, the phenomenon of the conical refraction was discovered for the elastic waves propagating 
T-H ■ along the acoustical axis of the crystal. For example, the inner conical refraction may take place when purely shearing 
waves propagate through the crystal along the third order symmetry axis (principal axis of symmetry [001] in triagonal 
'n ,; crystals; direction [111] along the diagonal of a cube in the cubical crystals). In certain cases phenomenon of the 
acoustical conical refraction may be rather pronounced. The angle of the conical refraction (i.e., angle of deflection of 
the beam from the acoustical axis of the crystal) may have significant magnitude, for instance, in calcit CaCO3«30°, 
Qh' in quartz ^17° etc. 

a' , Resuming aforesaid, the phenomenon of the conical refraction in the crystals (optical and acoustical) is known for 
quite a long time and hence is completely studied. However, further display of the phenomenon is rather surprising: 
^ ' in Ref. [1] there was theoretically discovered the phenomenon analogous to the inner conical refraction for magne- 
— H tosonic waves propagating in the isotropic plasma which is considered as a perfectly conducting fluid with adiabatical 
["T'" perturbations in it (see for brief description of the effect more accessible Ref. [2]). In particular [1], if the Alfven 
^ . velocity coincides with the speed of sound then in the case of propagating of a quasiplane magnetosonic wave along 
the uniform magnetic field the beams corresponding to the wave will form an empty cone with the opening angle 
53°08', i. e. the angle of conical refraction for the magnetosonic waves in isotropic plasma was predicted to be 26°34'. 
' In the present paper, while considering adiabatical disturbances in coUisionless plasma contained in a uniform exter- 
nal magnetic field and hence described by the theory of Chew, Goldberger and Low [3,4,5], is discovered phenomenon 
analogous to the inner conical refraction. Calculations showed that pending angle of the conical refraction for the 
magnetohydrodynamical (MHD) waves in the plasma with anisotropic thermal pressure should be 18°26'. The aim 
of this paper is to show theoretically that the phcmomenon analogous to the inner conical refraction should also exist 
in anisotropic plasma medium and to encourage other authors to prove experimentally existence of the phenomenon. 
The beams are interpreted as lines along which the group velocity vector is directed. 

In the concluding part of the paper there are discussed possibilities of experimental corroborating of the discovered 
phenomenon. 
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II. MAIN CONSIDERATION 



It is known [4,5] that the hnear MHD waves existing in the plasma described by the anisotropic thermal pressure 
obey to the following dispersion relation: 
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We further can rewrite eq. (1) in more convenient form introducing notations similar to those ones from Ref. [1]: 

a = r-. , c_L = W — , C||=W — 

where a is the Alfven velocity, c± and C|| are the sound speeds measured in perpendicular and parallel directions in 
respect to the uniform, external magnetic field Bq. Doing so, we get: 
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where V = iv/k denotes phase velocity of the waves. Traditionally, here denotes an angle between Bq and wave 
vector k. 

Let us consider the special subcase of eq. (2) assuming that condition c± ~ V^C|| is fulfilled in the medium under 
study. It is worthwhile to note that the latter assumption leads to existence of so called "angular point" [1] on the 
surfaces formed by the phase velocities (both fast and slow MHD waves) just in the point of intersection of the ones 
with the direction of the magnetic field Bq. Such a location of angular point, in turn, is responsible for the existence 
of the effect of conical refraction in the medim under study. The main guess was to find such a condition which would 
led to presence of the conical refraction in the anisotropic plasma medium, i.e. to discontinuous behavour of the 
vector of group velocity when tends to zero (see below). Since possessing such a condition we can rewrite eq.(2) as 
a following: 
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Basing on the latter expression we may calculate the group velocity of the waves. In general, the group velocity il 
has physical meaning of the wave packet propagation speed and can be found by simple formula u = dto/dk. The 
derivation of (3) results: 
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where we introduced unit vectors N = k/\k\ and A = Bo/\Bq\. To realize the physical picture geometrically it is 
worthwhile to state here following reasonings [1]: let multiply scalarly the latter formula by N: 



{uN) = V 
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i.e. projection of group velocity on the direction of wave nornal is equal to the phase velocity in the same direction. 
Multiplying eq.(5a) by k one can get: {uk) = u. Hence, for the infinitesimal changes of u, k and u the following 
relation is valid: {u5k) + (kSu) = 5u. According to the definition of group velocity {u5k) = 6u). Therefore, (kSu) = 0, 
or in other terms 
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The vector Su is located in the plane tangential to the surface (so called beam surface [1]) formed by the ends of 
vectors of the group velocities drawn away from the arbitrary point O. Therefore, eq.(5b) implies that the plane, 
tangential to the beam surface, is perpendicular to the corresponding wave nornal. i.e. the beam surface is the 
envelope for the plane wave fronts, which were propagated from reference point O in every directions per unit time. 
The latter statement in combination with eq.(5a) completely determines geometrical relationship between surface of 
phase velocities with the beam surface. However, it should be stressed that this geometrical relationship is valid only 
for the waves which phase velocities depend only on the direction and not on the length of the wave vector k. It is 
enough to show that in such case eq. (5a) is valid. Let V (phase velocity) be arbitrary function of wave vector k. 
Then u) = kV{k) and hence, u = dui/dk = V{k)k/k + kdV/dk. It is necessary and sufficient for condition (5a) to be 
hold that the V must satisfy the following requirment: kdV/dk = or that is the same as hidV/dhi = 0. Therefore, 
one can conclude that V must be homogeneous function of fcj components of the zeroth order. The latter statement 
is equivalent to the requirement that function V depends only on the direction of wave vector k and not on its length. 
As it can be seen clearly the dirpersion relations (1) and of course (2) and (3) satisfy this condition. 

It is known [1] that since possessing surface formed by phase velocities one can construct geometrically the beam 
surface (formed by the ends of vectors of the group velocities) and vice versa. But instead of doing this it more 
convenient to carry out analitical calculation. Let us describe a direction of the group velocity u by an angle (p 
between u and Bq (i.e., between u and A). One can determine the angle ip by the formula: 
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which with the aid of eq. (4) yields: 
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By means of the eq. (6) we completely determined the direction of the group velocity u whereas its module can be 
calculated by expression: 

\u\=V/cos{ip-e) (7) 

It is worthwhile to note that the eq. (7) can be obtained by using self-evident connection (5a) between the group 
velocity u and phase velocity V. 

It is simple to show that expressions (4) and (6) have uncertain behavior as 0/0 when a = c± (i.e., when a = c± = 
\/2c||) and 6 tends to zero. Removing the uncertainty yields: 

tan <p{+0) = lim tan (f{e)=±^ (8) 



u{+0) = (9) 

in eq. (8) sign "+" and "-" correspond to the slow and fast MHD waves respectively. The above treatment shows 

that when the wave normal N tends to reach a direction parallel to Bq (i. e. when 9 tends to zero), vector of group 
velocity u does not takes the same direction, but its limiting position (when 6* ^ 0) is the surface of an empty cone 
with opening angle 2(p{+0) = 2arctan(l/3) = 36°52'. When 9 is zero exactly then u is strictly parallel to Bq. Such 
a discontinuous behavior of vector u can be explained by means of existence of so called [1] angular point which in 
the case of the anisotropic plasma exists when a = c± = v^cij. (Note that for the waves in isotropic plasma the 
latter condition was rather simple [1] a = c where c denotes ordinary sound speed). So, when a = c± = \f2c\\ the 
angular point on the; surfaces formed by the phase velocities (both fast and slow MHD waves) is located in the point 
of intersection of the ones with the direction of the magnetic field Bq. Which, in turn leads to the discontinuous 
behaviur of u and causes the hole effect. In other terms the angular point is nothing else then point of intersection of 
the phase velocities of the fast and slow MHD waves. One can see from the Fig. 1 that when the condition a = c± 
is not satisfied than the polars of phase velocities (according to the dispersion relation (3)) do not have any points of 
intersection. It follows from the eq. (2) that when a = c± = \/2c|| , then 
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From the cq.(lO) we gather that when 9 = then V+ = V- and hence the phase polars have two points of intersection 
which can be clearly seen from the Fig. 2. It is apparent that in these points derivative of the phase velocities has 
uncertain behaviour. Removing of this uncertainity was done above. Result is clear now: in the anisotropic plasma 
for the linear MHD waves there should exist the effect analogous to the conical refraction and we worked out detailed 
explanation of the phenomenon. 

III. DISCUSSION 

Here in this section we discuss possibilities of experimental testing of existence of the phenomenon. Strictly speaking, 
the plane wave is idealization and real physical significance has quasiplane wave with small but finite angle of aperture. 
It should be emphasized that no real physical object corresponds to the wave with wave normal N exactly parallel 
to Bq [9 = 0), whereas in the case of propagation of quasiplane MHD wave with small but finite angle of aperture 
{6 0) the narrow beam of such quasiplane waves, as it was shown above, will form an empty cone with opening angle 
2(p{+0) = 36°52' if a = c_L = V^c\\ condition is fulfilled. We may suppose the following experiment for demonstration 
of the phenomenon analogous to the conical refraction: Let us consider two coUisionless plasma media contained in 
a uniform external magnetic field Bq. The media are separated by means of the plane perpendicular to the Bq. Let 
choose external conditions in the second medium so that to satisfy the condition a = cx = \/2c|| permanently. If 
we place a diaphragm between the first and second media then the beam of MHD waves propagating from the first 
medium towards the second one, just after passing the diaphragm will deflect from the direction of Bq and form an 
empty cone with angle of opening 2(p{+0) = 36°52'. 

It is worthwhile to note that the effect of the conical refraction may be used for plasma diagnostics. It is clear that 
the due to the effect of the conical refraction intensity of wave which is registered by a detector placed on the path of 
propagation of the wave will show steep drop for certain value of the magnetic field which satisfies the condition 

Such a technique was proposed for isotropic plasma case in Ref.[6] 

For the author is unknown whether the phenomenon analogous to the inner conical refraction of MHD waves in 
the anisotropic plasma was observed experimentally yet. However, one should take into account the difficulties in 
prepation of a plasma medium with permanently satified condition a = c± = a/2c||. The aim of this paper is to 
stimulate other authors, in spite of the techinacal difficulties, to prove the existence of the effect experimentally. 
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